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Ongoing interest in analysis and design of low sonic boom supersonic transports re- 
quires accurate and efficient Computational Fluid Dynamics (CFD) tools. Specialized grid 
generation techniques are employed to predict near-field acoustic signatures of these con- 
figurations. A fundamental examination of grid properties is performed including grid 
alignment with flow characteristics and element type. The issues affecting the robustness 
of cylindrical surface extrusion are illustrated. This study will compare three methods 
in the extrusion family of grid generation methods that produce grids aligned with the 
freestream Mach angle. These methods are applied to configurations from the First AIAA 
Sonic Boom Prediction Workshop. 


I. Introduction 

There is interest in developing business jet 1,2 and larger 3,4 classes of supersonic vehicles. Supersonic 
flight over land is considered crucial to the commercial success of these designs but is currently limited in 
many countries, including the United States of America. The acceptance of an aircraft’s sonic boom to the 
general population and regulatory agencies is a requirement for supersonic flights over land and therefore 
the commercial viability of a supersonic transport. Predicting how sonic boom signatures are perceived is 
a challenging task that requires the prediction of the signature on the ground. This task is complicated by 
long propagation distances, atmospheric variations, and the Earth’s turbulent boundary layer . 5 A detailed 
review of the history and state-of-the-art of sonic boom modeling is provided by Plotkin . 6 Surveys of sonic 
boom prediction methods are available from Ozcer' and Park . 8 

The propagation of a sonic boom is often separated into logical stages, or regions, to facilitate analysis . 6 
The near-held is a region near the aircraft where shocks are formed and are strongly influenced by nonlinear 
phenomena such as shock-shock interaction, shock curvature, and cross flow. Computational Fluid Dynamics 
(CFD) methods are typically employed in this near-field region. Sonic boom propagation methods 6 are 
typically employed beyond the near-field where the geometric details of the configuration are less important 
than atmospheric variations and molecular relaxation phenomena. These atmospheric propagation methods 
rely on the availability of an accurate near-field signature. 

Near-field CFD simulations often take advantage of the predictable characteristics and limited domains of 
influence of supersonic flow. There are many near-field CFD approaches and most of them include specialized 
techniques for grid generation and adaptation. A review of these previously published methods is provided 
in Section II, but there has been a limited number of comparisons between methods. The near-field CFD 
predictive capability of five methods was studied in an internal NASA workshop . 9 The American Institute 
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of Aeronautics and Astronautics (AIAA) is sponsoring the First AIAA Sonic Boom Prediction Workshop 
(SBPW-1) 10, 11 to examine methods employed by a broader section of the CFD community. 

Three methods in the family of grid extrusion methods are examined in this article. These methods 
extrude the outer boundary of a “core” grid, in a series of steps or layers, to produce grid elements that 
align to the freestream Mach angle /i, 

M = sin” 1 (1/M), (1) 

where M is Mach number. The extruded layers form a “collar” grid. These three methods will be applied 
to configurations from the SBPW-1 to document the properties of the methods. Sets of uniformly refined 
extruded grids were provided to the SBPW-1 for the two required cases. To introduce these extrusion 
methods, The effect of grid alignment on accuracy and the effect of extrusion surface topology on extrusion 
robustness are illustrated to indicate the key features of the extrusion process in Section III. 

II. Previous Work 

Cliff and Thomas 12 used structured and isotropic unstructured grid methods without alignment to the 
freestream Mach angle. Djomehri and Erickson 13 applied an adaptive isotropic unstructured grid FE- 
LISA 14, 15 method. A linear acoustics propagation code was required to reach the validation data locations 
for these isotropic methods. Fouladi 16 applied linear VGRID sources 17 and USM3D 18, 19 to isotropic grids 
of a body of revolution. Choi, Alonso, and van der Weide 20 applied isotropic grid adaptation in areas of 
large pressure gradient and produced good comparisons to distances less than one body length. Without 
alignment (using isotropic unstructured grids) only short distances were attempted or a propagation method 
was used to reach validation data locations. Carter and Deere 21 examined the effect of unstructured grid 
resolution and showed a significant improvement with alignment of anisotropic grid to the freestream Mach 
angle. These anisotropic alignment methods were improved and evaluated by Campbell et al. 22 Carter, 
Campbell, and Nayani 23 examined a low boom configuration with laminar and turbulent analysis. These 
freestream Mach alignment methods evolved into the extrusion methods of Nayani and Campbell 24 and Cliff 
et al. 25 

Meredith et al. 26 showed an overset grid system with Cartesian and freestream shock-aligned grids for the 
Shaped Sonic Boom Demonstrator (SSBD) vehicle. The shock-aligned grid provided better comparison to 
flight data measurements. Bui 27 studied the grid resolution and alignment. He found that the Mach-aligned 
grid provided better prediction of cone-cylinder measurements than a Cartesian aligned grid, even at much 
higher spacial resolutions. Ishikawa et al. 28 utilized structured grids inclined to Mach angle and examined 
the effect of grid resolution. 

Siclari and Darden 29 detailed efforts to align structured grids and adapt to shock locations. This adaptive 
approach was also employed by Page and Plotkin. 30 Kandil et al. 31,32 used various Euler methods very near 
the body. This near body solution was interpolated to an adaptively shock-aligned structured grid. Complete 
details of the shock-alignment method is provided by Ozcer ' 

Ozcer and Kandil 33 applied FUN3D 34,35 with an OptiGRID 36 anisotropic grid adaptation method to 
resolve the pressure gradient by refining and aligning the tetrahedral grid. Alauzet and Loseille 3 ' utilized 
anisotropic metric-based interpolation control of Mach number. The anisotropic adaptive method allowed 
second-order global mesh convergence in the Ly-norm of Mach field interpolation error. This framework 
was extended to goal-oriented adaptation by Loseille, Dervieux, and Alauzet. 38 Jones, Nielsen, and Park 39 
adapted grids to a Mach number Hessian grid metric scaled with an off-body pressure integral output-based 
grid resolution request. The use of tetrahedral grid allowed arbitrary alignment to local flow characteristics. 
The output-based adaptation scheme was extended to tetrahedral cut-cells 8 and viscous effects with frozen 
boundary layer grids. 40 Barter 41 applied an output-adaptive higher-order discontinuous Galerkin scheme 
with artificial viscosity to a laminar airfoil. Fidkowski and Darmofal 42 developed a cut-cell scheme and 
adapted to an off-body pressure integral. 

Nemec, Aftosmis, and Wintzer 43 used Cart3D 44 along with adjoint-based mesh adaptation to predict 
the pressure signature of a diamond airfoil. Improvements due to grid alignment to freestream Mach angle, 
cost function formulation, and aspect ratio were examined. Cliff et al. 45 illustrated the benefit of aligning 
a Cartesian direction with the freestream Mach angle. Aftosmis et al. 46 showed that pressure signatures on 
the adapted grids at different off-track angles collapse to the same signature, as expected in axisymmetric 
flow. 

Adaptation is important for complex flow, but a priori grid alignment has produced satisfactory results 
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for problems with know characteristics. Alonso, Kroo, and Jameson 4 ' generated structured H-mesh with 
grid lines canted at the freestream Mach angle to maximize the resolution of the pressure signature. The grid 
alignment has also been implemented by interpolating between multiple solvers or grid constructions. Laflin, 
Klausmeyer, and Chaffin 48 used an adapted isotropic grid very near the model to facilitate grid generation. 
The NSU3D 49 solver computed the solution in this region, which was interpolated to a structured mesh that 
is swept to align mesh cell faces with the freestream Mach lines and solved with CFL3D. 50 

Waithe 51 compared USM3D 49 with generated isotropic and adapted anisotropic unstructured grids to 
OVERFLOW 52 with a freestream shock aligned overset grid. A hybrid method was employed that inter- 
polates the unstructured solution to a cylinder for OVERFLOW initialization. Howe 53 used Cart3D in 
combination with an OVERFLOW freestream shock aligned overset grid and examined how the location 
of interpolation surface affects the predicted signature. Ishikawa et al. 54 utilized a hybrid method with an 
inner unstructured grid solution interpolated to a Mach angle aligned structured grid. This hybrid approach 
was also utilized by Morgenstern, Buonanno, and Marconi. 55 


III. Grid Generation Considerations 

This article will extend previous efforts 21, 27 ’ 43 > 4 5 examine the alignment of the grid to the flow features. 
Specifically, the previous works compared Cartesian aligned or isotropic unstructured grids to freestream 
Mach aligned grids. The current work investigates the degradation due to a small degree of misalignment 
that can result for complex flows around realistic configurations. 

The issues related to extruding a cylindrical boundary introduced by Nayani and Campbell 24 and Cliff 
et al. 25 is further investigated in this article. These methods extrude the outer boundary of a core grid to 
create a collar grid. Section III.C details the conditions that can result in inverted prismatic elements that 
invalidate the extruded collar grid and strategies to address this pitfall. Geometric growth is used by these 
methods to gradually increase the height of subsequent layers of this extraction. 


III. A. Geometric Growth 

A geometric growth of layer thickness is used by the extrusion methods. The first layer is typically constructed 
to approximately match cells sizes at the inner core grid and collar grid interface. The number of elements is 
reduced by increasing layer thickness and cell aspect ratio away from the model. Geometric growth produces 
a smooth variation of cell heights from the first layer to the last layer to accomplish these goals. The initial 
prism altitude or thickness to and growth rate r determines the thickness of layer i, 

U = t 0 r\ (2) 


where i is 0, . . . , n — 1 and n is the total number of layers. The total thickness T is given by, 


and Newton’s method is used to solve for r when t 0 , n, and T are specified. 


( 3 ) 


III.B. Grid Alignment 


A hexahedral grid consisting of a single span of cells is utilized to examine the impact of grid alignment on 
propagated shocks. This yields a 2D solution. The lower boundary is given by the expression, 


s(z) = < 


0 

O.OLr 

0.01 

O.OLr - 0.01 
0.02 


x < 0 

0 <= x <= 1 

1 <= x <= 2 

2 <= x <= 3 
0 < x 


( 4 ) 


These two ramps produce two pairs of oblique shocks and expansions as shown in Fig. 1. This combination 
of shocks and expansions is expected to transform from square wave into a sawtooth wave as the signal 
propagates, as seen in Fig. 2(a). The grid is extruded at the freestream Mach angle in Fig. 2(b) and has a 
constant spacing in the r-direction of h x = 0.125. The grid has a first cell height of to = 0.1 and growth 
rate of r = 1.051. The total height is T = 100 after 80 layers. 
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Figure 1. 2D solution diagram. 




(a) Sheared domain and pressure disturbance. (b) Close-up of grid sheared at freestream Mach angle. 


Figure 2. 2D sheared grids. 


III.B.l. Signature Propagation for Aligned and Misaligned Grids 

Two additional grids are created with the shear angle four degrees ahead and four degrees behind the 
freestream 1.6 Mach angle. The signatures produced by the three grids sheared ahead, at, and behind the 
freestream Mach angle are studied to observe the impact of grid alignment. The flux- vector splitting of 
van Leer 56 is used with unlimited unweighted least-squares reconstruction as described by Anderson and 
Bonhaus. 34 The pressure sampled at increasing propagation distance is shown in Fig. 3 for the three grids 
sheared ahead, at, and behind the freestream Mach angle. For H = 16.0 ( H is the displacement in the 
^-direction) and less the three grids produce very similar signatures as the pressure disturbance transforms 
from square waves to sawtooth waves, see Fig. 3(a-c). At H = 32.0, a difference in the three shear angles is 
observed, see Fig. 3(d). The Mach angle aligned grid is producing the expected sawtooth signature and the 
misaligned grids are producing a smoothed signature. The signature of the grid that is inclined four degrees 
ahead of the Mach angle is more smooth than the signature of the grid that is inclined four degrees behind 
the Mach angle. 

III.B.2. Element Type and Refinement of Aligned and Misaligned Grids 

The effect of grid refinement and element type for Mach angle aligned and nonaligned grids at H = 64 
is shown in Fig. 4. The aligned grid exhibits the expected sawtooth signature for all configurations. The 
smooth signatures at H = 32.0 of nonaligned grids is further smoothed at H = 64. The two sawtooth 
waves have merged into a single smooth wave for the four degrees ahead of Mach angle grid, see Fig. 4(a). 
Converting the hex grid into a single layer of prisms also produces a single smooth signature for the four 
degrees behind of Mach angle grid in Fig. 4(b). Grid refinement is performed in the x and 2 directions. The 
geometric growth rate for the refined 3 grid is r = 1.025 to yield 160 layers with a total height of T = 100. 
The x refinement, shown in Fig. 4(c), does not significantly improve the sharpness of the nonaligned grids, 
but z refinement, shown in Fig. 4(d), does improve the nonaligned grid signatures. The improvement with 
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dp/pinf dp/pinf 


4 Degrees Ahead A 4 Degrees Ahead 




(a) Pressure at H = 0.05. 


(b) Pressure at H = 8.0. 


4 Degrees Ahead 



4 Degrees Ahead 



(c) Pressure at H = 16.0. 


(d) Pressure at H = 32.0. 


Figure 3. 2D sheared grids. 
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z refinement may be connected to a reduction in cell aspect ratio for z refinement, which appears to benefit 
the misaligned grids. 


4 Degrees Ahead 



4 Degrees Ahead 



(a) Original spacing, hexes. 


(b) Original spacing, prisms. 


4 Degrees Ahead 



4 Degrees Ahead 



(c) Refinement in the fc-direction, 0.5 h x - 


(d) Refinement in the ^-direction, 0.5 h z . 


Figure 4. Pressure at H — 64.0. 


III.C. Cylindrical Surface Extrusion Pitfalls 

The three methods described in this article all extrude elements aligned with the freestream Mach angle 
from a triangulated surface. Even when the the nodes of this triangulation are projected to a cylinder, the 
normals of the linear triangles deviate from the analytical normal of the cylinder. If the dot-product of the 
triangle normal and extrusion direction becomes negative, an inverted prism is created that invalidates the 
grid. 

Shewchuck 57 examines the approximation of paraboloid surface normals by linear triangles. He illus- 
trates that a triangulation with large angles have alternating normal directions like a “washboard.” Poor 
approximation of a cylindrical surface produces surface normals that are not normal to the cylinder axis. 
These inclined triangles are used to form the bases of prism volumes that are extruded in the direction 
of freestream Mach angle. The radius of the cylindrical surface grows with successive extrusions and the 
triangles are elongated in the tangential direction. The elongation of a triangle with each node on a cylinder 
at different tangential locations will result in a decrease of the angle between the triangle’s normal and the 
a cylinder axis. When the triangle’s normal is inclined more than the freestream Mach angle, an extruded 
prism may invert. 

Further study into the cause of the negative cells has shown that this is an inherent feature for grid 
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extrusion methods if the core outer grid has surface triangles with a downstream point that does not line up 
directly behind one of the upstream points. The special case for a triangle with one side perpendicular to 
the cylinder axis is illustrated in Figure 5. All 3 points of the triangle are on the outer cylinder of radius R 
and point c is downstream of points a and b by a distance dx. Points a and b are offset from point c by a 
radial angle </). This offset results in a normal distance dz which is given by 

dz = i?(l — cos ((/>)) (5) 

resulting in a face inclination angle A computed as 

A = taxi~ 1 (dz/dx) (6) 


X 



Front View Side View 


Figure 5. Illustration of non-aligned outer boundary face inclination problem. 

As layers of cells are added to the collar grid, dx remains fixed, but the value of R increases, thus 
increasing dy and A. This is shown by the triangle with points A, B, and C. The A, B, and C triangle is 
only twice the radius of the a, b, and c triangle and A has increased significantly in the “Side View.” The 
R can increase twentyfold in a typical application, which would increase A by a greater amount. Eventually, 
the boundary face inclination angle becomes parallel to the interior cell faces that are sheared at the Mach 
angle and the volume of the cell connected to that face becomes negative. 

There are methods to control the angle between the surface normal and the cylinder axis. Arbitrary 
triangles with small angles will reduce this deviation, 57 but not eliminate it. Constraining each extruded 
surface triangle so that one side is parallel to the cylinder axis will eliminate the error in the x-r plane. A 
practical method to produce both triangles with an edge aligned to cylinder axis and maintain small angles 
is subdividing a structured grid, aligned to the curvature direction, into right triangles. A structured grid 
aligned to the curvature constrains variations in grid density. To allow a variation in grid density, a prism 
with a regular polygonal base can replace cylindrical surface. This creates an approximation to the cylinder 
with planar surfaces. The amount of approximation can be controlled by the number of polygonal sides. 
If the topology of this ruled surface between the polygonal bases is enforced during grid generation, the 
normal of planar surfaces is expressed exactly by any triangle regardless of shape and each triangle is exactly 
normal to the center axis by construction. Without constraints on the construction of the extrusion surface, 
refinement that splits long edges and can improve grids, but does not guarantee robustness. 
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IV. Extrusion Methods 


Three methods are introduced to extrude the outer boundary of a core grid into layers that form a collar 
grid aligned to the freestream Mach angle. The geometric growth, grid alignment, and extrusion robustness 
discussions in the previous section are used to aid the description of these methods. 

IV. A. Boom Grid (BG) Generation Method 

The Boom Grid (BG) method 24 is similar to other “extrusion” grid generation approaches in that the new 
grid is created by extending lines from grid points on a boundary surface of an existing grid in a direction 
approximately normal to the surface. The new cells are extruded through the outer, typically cylindrical, 
farfield boundary of an existing “core” 3D grid around a configuration, creating a “collar” grid, see Fig. 6. 
For triangular outer boundary faces on the core grid, the extruded cells are prisms. For methods that only 
support tetrahedra, each of these prisms is split into three tetrahedral cells. The BG method utilizes very 
efficient cell generation and splitting algorithms that require a few seconds to generate hundreds of millions 
of new tetrahedral collar grid cells. 



Figure 6. Prism extrusion process in the BG grid methodology. 


The extension direction extends radially from the reference point when viewed from the front, see Fig. 6. 
The extension direction is swept to the Mach angle when viewed from the side. This grid-line shearing, 
along with stretching of the distance between subsequent layers, creates high-aspect ratio cell faces that are 
closely aligned with the Mach angle so that the numerical dissipation of the boom signature is reduced. As 
few as 20 or more than 100 layers are generated to reach the desired location for extracting a sonic boom 
signature. 

If the outer boundary is a cylinder and the collar grid lines are extruded normal to the original core grid 
outer boundary, then all of the extruded grid lines appear to originate at the center of the circle (front view 
of the cylinder) as a radial reference point. If the core grid has the configuration located at this reference 
point, then boom signature propagation lines will be approximately aligned with the radial grid lines, as 
desired. The weakness of this approach, however, is that there can be a significant distance from the bottom 
of the fuselage to the beginning of the collar grid. As the grid is not sheared and stretched in this core grid 
region, the signal can be dissipated before reaching the collar grid. 

If only on-track signatures are needed, then the above problem can be resolved (though with some 
addition cost due to more grid cells) by simply increasing the radius of the outer core grid cylinder to allow 
the configuration to be shifted down closer to the lower boundary. For off-track signatures, however, this 
approach gives poor alignment of the grid lines and signal propagation direction. The BG method attempts 
to address this problem by using a radial reference point that is not the center of the circle. This allows the 
configuration to be shifted down closer to the outer boundary and still have the radial lines that emanate 
from the configuration for better signal propagation. The radial reference point can be located to focus the 
grid lines in a desired off-track direction to provide better signal resolution at the off-track angle with very 
little degradation of the on-track signature. Figure 7 illustrates this capability, showing collar grid inflow 
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planes for standard on-track and focused off-track cases. Nayani and Campbell 24 provides more detail on 
this option. 



(a) On-track. 


(b) Off-track. 


Figure 7. 


Effect of BG grid-focusing option for off-track analysis on collar grid inflow plane. 


As noted by Nayani and Campbell, 24 the BG method is an improvement over the previous SSGRID 22,58 
and SSG 24 grid stretching and shearing codes in terms of being able to reach a given distance below the 
aircraft before creating negative volume cells. BG has successfully generated collar grids with over 1,400 
cell layers to distances of 100 body lengths without negative cells when the original core grid met the 
requirement of every face having one edge parallel to the cylindrical outer boundary axis. An automated 
method for creating an input file for the VGRID grid generation method that enforces this requirement has 
been developed and is referred to as the Cylindrical Outer Boundary (COB) code. In addition to achieving 
large signature location distances, the use of COB also allows the BG code to successfully obtain grids 
focused at off-track angles up to 90 degrees. 

IV. B. Mach Cone Aligned Prism (MCAP) Generation Method 

The Mach Cone Aligned Prism (MCAP) collar grid approach was developed by Cliff et al. 25 and applied 
the to Delta Wing Body and Low-Boom Wing Body Tail configurations. Elmiligui et al. 59 applied MCAP 
to a wing body configuration. MCAP extrudes the triangular faces of the cylindrical outer boundaries of 
a core tetrahedral mesh into a prismatic collar grid where the prisms are subdivided into tetrahedra. The 
shearing angle of the prism cells is adjusted for the angle of attack and Mach angle to allow for both on- and 
off-track aligned grids. These prisms are subdivided into tetrahedra to produce a grid composed entirely 
of tetrahedral cells. This provides an automated process to construct grids suited to obtain accurate sonic 
boom pressure signatures. Thus, configurations with any level of geometric complexity can be analyzed, 
provided a tetrahedral near-body volume grid of tetrahedra with a cylindrical boundary is supplied. 

The core near-body grids may be constructed in a manner that is suitable for either viscous or inviscid 
computations. The radius of the cylinder should be small so that the boundaries are close enough to the 
vehicle to allow the Mach cone aligned prism cells to begin near the body before signals are dissipated by 
the near-body tetrahedra. An example of a cylindrically-shaped mesh and configuration is shown in Fig. 8. 
The fore and aft boundary faces and symmetry plane are removed in the figure for clarity. These fore and 
aft boundaries can be flat or aligned with the Mach angle. Shearing the aft boundary to form an inverted 
cone shape is an effective technique to remove unnecessary points downstream of the model. 

The mesh density within the cylindrical grid needs to be sufficiently dense within the sonic boom region 
of influence so that the pressure signature can be taken far from the model, and configuration variations in a 
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Figure 8. Cylindrical Mesh. 


design environment would not require different volume grid densities. The grid will need to be highly refined 
within the Mach cone and behind the model for approximately a third to a half body length to capture 
details of the bow, tail, and return to ambient pressure signatures. Typically, refinement is necessary up 
to 90 degrees off-track even for on-track computations. To produce a high-density collar grid, a refined 
triangular cylinder grid is required. 

MCAP inherits the grid spacing in the axial direction from the cylinder triangulation, and adjusts the 
radial stretching and shearing to align with the Mach cone angle around the aircraft. Each triangular face 
is extruded to form a prism that is sheared to align with the Mach angle. Each prism is split into three 
tetrahedra, which adds an edge to each of the quadrilateral faces of the prism cells. The quadrilateral face 
of the prism may not be a parallelogram since the extruded triangle edge is larger in the tangential direction 
than the previous layer due to the increasing radius of the cylinder. A prism and the three tetrahedra from 
splitting the prism are shown in Fig. 9. The topology of the subdivision of each prism into three tetrahedra 
is the same for all layers, as shown in a shaded and exploded view in Fig. 9(b). 




(a) Extruded prism. (b) Subdivision into tetrahedra (exploded view). 

Figure 9. Prism extrusion and subdivision. 


The extruded cylindrical mesh is sheared downstream at each layer for alignment of the mesh with the 
freestream Mach cone at angle /x. Each mesh point is shifted by 


dx 


t 

tan(ix) ’ 


( 7 ) 


where dx is added to the x coordinate shown in Fig. 10 and t is the thickness of the layer in the radial 
direction. When the angle of attack is zero, v = /x, which results in a symmetric shear. The thickness for 
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each successive layer can grow geometrically, as shown in Fig. 10, or be held constant if the growth ratio is 
set to one. 



Figure 10. Core and extruded collar symmetry plane grid and model. 


When the angle of attack a is not zero, then an additional asymmetrical shear is applied that depends 
on a and (f>, where <f> = tan _1 (— z/y). A periodic weighting of a is given by 

( 8 ) 

and the shear angle becomes 

v = n + (— a(l — a) + aa). (9) 

The weighting function is an approximation to the intersection of an a-rotated cone and a cylinder where 
the intersection is an ellipse in the a-rotated plane. This approximation is sufficient for supersonic vehicles, 
where angles of attack are expected to be small. 

With each additional prism layer there is increased stretching of the mesh in the tangential direction 
due to the growth of the radius to the outer boundary. Thus, in order to obtain accurate solutions at 
large distances from an aircraft the stretched triangles and tetrahedra should be split to improve tangential 
refinement of the mesh. Mesh splitting should only be used when necessary because it approximately doubles 
the mesh point count and triples the cell count on the split layer and all successive prism layers. 

The highly stretched tetrahedra on the outer prism layers have faces with two very long edges and one 
short edge. Thus, in order to reduce the tangential stretching of the prism mesh, the MCAP software allows 
the user to split the tetrahedra of one or more prism layers during the construction of the prismatic mesh. 
The algorithm splits only the longest two edges of each face on the outer boundary and the tetrahedra 
associated with this refined face. Splitting all edges would have been less difficult to program but would not 
have alleviated the stretching or improved the quality of the mesh. 

Splitting amounts to adding data to lists of points, faces, and tetrahedra. The algorithm loops over 
the indices of the face elements on the outer cylinder. Each face is typically split into three faces with the 
same connectivity as shown in Fig. 11. The labeled points as described in the figure text can be found 
regardless of orientation. The algorithm determines whether the edges of neighboring faces had been split 
before splitting. Therefore if a small edge of a previous split triangle is common with the mid-length or long 
edge of the current face, then that edge will not be split. Thus it is possible to not split the face or split it 
into 2 to 4 triangles. There are a finite number of conditions encountered and the connectivity of face points 
is defined depending on the condition. 

The tetrahedra associated with the split faces must be split as well. This involves two searches of the 
tetrahedron list. First, search for the tetrahedron with three common points of the face being split. The new 
tetrahedra are simply the faces plus the fourth point of the un-split tetrahedron illustrated by dashed edges 
in Fig. 11. The second search is to find any tetrahedra that contain the fourth point and two points of the 
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PI Common point between largest and mid length edges 

MP1 Midpoint of longest edge 

P2 Point opposite mid length edge 

MP3 Midpoint of smallest edge 

P3 Point opposite longest edge 

MP2 Midpoint of mid length edge 



Figure 11. Details of extrusion layer refinement. 


face being split, shown by the center medium gray cell in the example. If the edge of the abutting triangle 
is split then the neighboring tetrahedron must be split. This is not the case in this example, but would 
be if the upper dark gray tetrahedral were split at MP3. These two searches are typically non-exhaustive 
since they terminate when the tetrahedron is found. Each face and associated tetrahedron are split until 
the face list on the outer cylinder is exhausted. The program then copies the refined layer connectivity for 
the remaining layers of prismatic cells. Each newly added point is adjusted radially so that it is the same 
distance from the cylinder axis as its neighbors. This adjustment reduces the deviation between the triangle 
normal and the normal of a cylindrical surface by reducing dy in Fig. 5. 


IV. C. Inflate Generation Method 


The Inflate collar grid generation method can utilize a core grid with a triangular surface from any source. 
To guarantee robustness, the outer core grid surface should be composed of planar surfaces parallel to the 
cylinder axis. A prism with regular polygonal bases satisfies this constraint and the sides of the regular 
polygons can be increased to improve the cylindrical approximation of the ruled surfaces. 

To form a layer, each node in the current triangular surface is first extruded in the radial direction d 
a distance l so that the dot product of the displacement vector and the adjacent triangle normals is the 
requested thickness fj, 

l = dr • n, (10) 


d l dr + tan(iz) 


dx. 


( 11 ) 


where n a unit triangle normal, dr is a unit vector in the radial direction, and dx is a unit vector in 
the positive x direction. Calculating the extrusion distance of each node maintains planar surfaces in each 
successive layer. This is accomplished through Eq. (11) and Eq. (10) where the extrusion distance is increased 
as the adjacent triangle normals deviate from the analytic normal of a cylindrical surface. 

An example of a topology with regular polygonal outer boundary is shown in Fig. 12(a). It is a quarter 
domain of with a cone cylinder in the lower right portion of Fig. 12(a). Symmetry planes for the bottom 
and right of Fig. 12(a). The outer boundary in the upper left is composed of ruled surfaces connecting the 
sections of regular polygons at either end. In this example, four sides are shown of the twelve sides that 
compose the entire polygon. The outer boundary is a poor approximation of a cylinder in this example, but 
the number of sides can be arbitrarily increased to reduce the approximation. Each of these ruled surfaces 
is extruded to create the topology in Fig. 12(b). 
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(a) Core grid topology. (b) Core and inflated collar grid topology. 

Figure 12. Topology of a quarter of a cone-cylinder model domain. 


The surface grids of the typologies shown in Fig. 12 are shown in Fig. 13. The core grid in Fig. 13(a) 
is composed of tetrahedra and the outer surfaces are composed of triangles with arbitrarily resolution and 
orientation. Each outer surface triangle in Fig. 13(a) is extruded into a prism in the same fashion as Fig. 9(a). 
The quadrilateral sides of the prisms are seen in Fig. 13(b) and each layer remains parallel to the original 
core grid outer surfaces. 



(a) Core surface grid. (b) Core and collar inflated surface grid. 

Figure 13. Surface grid of a quarter of a cone-cylinder model domain. 


V. First AIAA Sonic Boom Prediction Workshop Grids 

A series of uniformly refined mixed-element and purely-tetrahedral unstructured grids is developed for 
SBPW-1 with the following considerations: It is expected that the participants would use both cell-centered 
and node-centered finite volume methods, which have different resolution requirements. 60 The AIAA Drag 
Prediction Workshop series emphasized the benefit of a grid convergence study in interpreting results. 61 
Some methods require purely-tetrahedral grids and other methods prefer mixed-element grids. Determining 
impact of grid element types is also of interest. These considerations are addressed by a series of grids with 
consistent spacing distributions and common node locations between mixed-element and purely-tetrahedral 
grids. The Inflate method is used to create the required unstructured mixed-element and purely-tetrahedral 
grids for SBPW-1. 
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V.A. SEEB-ALR Model 


The SEEB-ALR model is described by Morgenstern et al. 4 It is designed as an axisymmetric body, but a 
very slight asymmetry was introduced during wind tunnel model construction. The SEEB-ALR solid model 
is constructed in NX® from a STEP file of the scanned as-built wind tunnel model. A 36-sided prism 
with regular polygonal bases is created by constructive solid geometry via a batch mode of the GridEx 
framework. 62 This prism is as small as possible while containing the SEEB-ALR model to minimize the use 
of isotropic grids that are not aligned to flow characteristics. The solid model of the SEEB-ALR is Boolean 
subtracted from the constructed prism. The triangular surface and tetrahedral volume grid are constructed 
in GridEx using the VGRID size specification 17 and the FELISA mesher. 14,15 

The surface grid after inflation is shown in Fig. 14(a). The rows of quadrilaterals in Fig. 14 are the sides 
of prisms that result by extruding the outer triangular boundary of the core grid into prisms. This is the 
coarsest grid h = 2.0 in the uniformly refined series and there is only one prism cell for each planar surface in 
Fig. 14(a). The configuration symmetry plane is shown in Fig. 14(b). The extrusion direction is specified by 
the Mach number to align to the freestream Mach angle and the layer thickness grows geometrically, Eq. (2). 
The thickness of the first extruded layer is intended to approximately match the size of adjacent tetrahedra. 
The core grid is isotropic with fine resolution at the nose, see Fig. 14(b). 



(a) Surface grid for the inflated grid, h = 2.0. (b) Inflated grid. 


Figure 14. SEEB-ALR model with inflated patch topology and surface grid. 


A family of uniformly refined grids was developed for the workshop and these grids are required cases. 10 
The total thickness of the inflation is T = 60 inches. The VGRID spacing field of the inner grid is multiplied 
by h and a series of grids are produced as listed in Table 1 with the inflation parameters. The r is computed 
to have a T = 60 as the number of layers n is varied by inverting Eq. (3). The complete generated grid size 
of the mixed-element grid is listed in Table 2. The prisms of the mixed-element grids are subdivided into 
tetrahedra 63 to form the purely-tetrahedral grids, Table 3. 


Table 1. SEEB-ALR inflation parameters. 


Characteristic h 

*0 

r 

n 

0.80 

0.04 

1.031 

125 

1.00 

0.05 

1.040 

100 

1.25 

0.0625 

1.050 

80 

1.56 

0.078 

1.063 

64 

2.00 

0.1 

1.081 

50 
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Table 2. SEEB-ALR mixed-element grids. 


Characteristic h 

Nodes 

Triangles 

Quadrilaterals 

Tetrahedra 

Prisms 

0.80 

10,159,421 

430,082 

193,875 

12,015,051 

15,850,125 

1.00 

5,252,466 

276,490 

123,100 

6,164,094 

8,161,100 

1.25 

2,537,977 

172,098 

80,080 

3,226,983 

3,822,800 

1.56 

1,165,851 

106,222 

50,240 

1,614,796 

1,682,496 

2.00 

674,734 

68,786 

31,250 

770,199 

1,023,350 


Table 3. SEEB-ALR, purely-tetrahedral grids. 

Characteristic h Nodes Triangles Tetrahedra 

0.80 

10,159,421 

817,832 

59,565,426 

1.00 

5,252,466 

522,690 

30,647,394 

1.25 

2,537,977 

332,258 

14,695,383 

1.56 

1,165,851 

206,702 

6,662,284 

2.00 

674,734 

131,286 

3,840,249 


V.B. Delta Wing Body Model 

Originally described as Model 4 in Hunton, Hicks, and Mendoza, 64 the solid model of the Delta Wing Body 
is extracted from a STEP file exported from CATIA® V5, where it was created from the description in the 
original report. The sting was created from drawings used to machine the wind tunnel sting. A simplification 
was made by omitting the strain gauges, wires, and instrumentation slots visible in a photograph of the model, 
see Fig. 15. The sting was inserted 0.079184 meters (3.1175 inches) into the wing body and a Boolean addition 
produced the solid model. 

The core region extends well beyond the fuselage to encompass the Delta Wing Body wing, see Fig. 17. A 
variation in triangle size can be in Fig. 17(b) for the h = 2.0 outer core grid surface. The leading edge of the 
extrusion planes are three prisms wide for most planes, see Fig. 18(b). The symmetry plane grid is shown in 
Fig. 19. The core grid is isotropic and randomly oriented, see Fig. 19(a). The initial layer thickness is twice 
the adjacent core grid spacing, see Fig. 19(b). The grid parameters are listed in Table 4. The portion of the 
mixed-element grid that is prismatic is smaller for the Delta Wing Body than the SEEB-ALR, see Table 2 
and Table 5. 



Figure 15. Delta Wing Body and sting. 


VI. Signatures 

A flow solution is computed on the grids generated by the three methods and pressure signatures are 
extracted. USM3D 18,19 is applied to the BG and MCAP grids. The FUN3D 34,35 flow solver is applied to 
the Inflate grids; the van Leer 56 flux is used with van Albada limiter 65 for this Euler solution. Examining 
the signatures from the uniformly refined mixed-element Inflate grids are intended to demonstrate grid 
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Figure 16. Delta Wing Body surface, h — 2.0. 




(a) Delta Wing Body model and polygon topology. 


(b) Regular polygon triangular surface, h = 


Figure 17. Delta Wing Body model with regular polygon patch topology and surface grid. 


Table 4. Delta Wing Body inflation parameters. 


Characteristic h 

*0 

r 

n 

1.00 

0.002500 

1.022 

100 

1.25 

0.003125 

1.028 

80 

1.56 

0.003900 

1.035 

64 

2.00 

0.005000 

1.045 

50 
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(a) Inner grid. (b) Inflated grid. 

Figure 19. Delta Wing Body model surface grid with the symmetry plane of the original and inflated grids. 


Characteristic h 

Table 5. Delta Wing Body mixed-element grids. 

Nodes Triangles Quadrilaterals Tetrahedra 

Prisms 

1.00 

5,363,065 

249,640 

55,000 

17,734,283 

4,689,000 

1.25 

2,742,341 

159,652 

34,880 

9,081,579 

2,369,920 

1.56 

1,411,515 

102,688 

22,464 

4,671,872 

1,205,440 

2.00 

695,709 

63,226 

13,700 

2,217,575 

613,000 
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Table 6. Delta Wing Body purely-tetrahedral grids. 


Characteristic h 

Nodes 

Triangles 

Tetrahedra 

1.00 

5,363,065 

359,640 

31,801,283 

1.25 

2,742,341 

229,412 

16,191,339 

1.56 

1,411,515 

147,616 

8,288,192 

2.00 

695,709 

90,626 

4,056,575 


convergence. The finest mixed-element Inflate, finest purely-tetrahedra Inflate, purely-tetrahedra BG and 
purely-tetrahedra MCAP are combined with experimental data to illustrate the properties of each method. 

VI. A. SEEB-ALR Signature 

FUN3D 34 ’ 35 is applied uniformly refined mixed-element Inflate grids for the SEEB-ALR. The signature at 
H = 21.2 inches is shown in Fig. 20. All five mixed-element grids are shown in Fig. 20(a). The similarity for 
the signatures from these series of refined grids indicates that grids are sufficiently resolved. Only a slight 
increase in sharpness at the extrema is observed. 




(a) Uniformly refined mixed-element grids (b) Various grid methods. 

Figure 20. SEEB-ALR, signatures at H = 21.2 inches 


The MCAP and BG grid signatures are plotted with the finest h = 0.80 mixed-element and the finest 
purely-tetrahedral grid signatures in Fig. 20(b). The purely-tetrahedral and mixed-element Inflate grid sig- 
natures over plot on this scale. The MCAP-Tet-USM3D and BG-Tet-USM3D also produce similar signatures 
with a higher pressure flat-top and lower pressure expansion than the Inflate grids. The reason for this slight 
difference is unknown. The model surface grids are inspected to eliminate geometry errors as a source of 
the difference. The amount of difference between the signatures is placed in context by comparing it to the 
magnitude of the experimental uncertainty. 

The sonic boom pressure signatures were measured in the NASA Ames 9- by 7-Foot Supersonic Wind 
Tunnel at Mach 1.6 and 1.7 in 2011 and 2012 during phase I and II of the NASA High Speed Project 
N+2 Experimental Systems Validations NASA Research Announcement (NRA) contract. 4 The sonic boom 
test data were obtained using pressure rail instrumentation developed by NASA in 20 1 0. 66 The rail design 
effectively eliminates amplification or reflection of the shock wave strength at pressure measuring orifices on 
the tip of rail, and is designated the reflection factor 1.0 (RF1) rail. Spatially averaged pressure signatures 
measured at multiple models positions along the RF1 rail were used to provide the accurate sonic test data. 67 
Additional description of the experimental data for a variety of flow conditions, including the data near the 
flow conditions of the SBPW-1 test cases, are provided in a companion paper. 68 

The computed signatures are compared with spatially-averaged experimental measurements from two 
sets of wind tunnel test runs in Fig. 21. All methods are within the experimental uncertainty for the flat-top 
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portion of the forward signature x = [28,37] and below the experimental expansion pressure x = [38,43] 
including uncertainty. The final recovery pressure x = [45, 55] is within the uncertainty of the experimental 
measurements. 



Figure 21. SEEB-ALR signatures at H = 21.2 inches, various grid methods and experimental data. 68 


VI. B. Delta Wing Body Signature 

The SEEB-ALR is an axisymmetric body (with a small eccentricity due to as-built tolerances). The Delta 
Wing Body has axisymmetric body with a wing, so it provides axisymmetric signatures forward of the wing 
and with varying off-track signatures aft of the wing for evaluating the suitability of grids for off-track 
propagation. The signature at H = 0.5 inches below the model is shown in Fig. 22. All four mixed-element 
grids are shown in Fig. 22(a). The signature changes little with grid refinement in the forward portion of the 
signature. The model base region x = [0.18, 0.28] changes the most with grid refinement and the faceting of 
the sting surface grid is evident for the coarser grids aft of x = 0.28. 




(a) Uniformly refined mixed-element grids. 


(b) Various grid methods. 


Figure 22. Delta Wing Body signatures at H = 0.5 and 0 = 0. 


The BG grid has clustering in the <f> = 0 direction as is shown in Fig. 7(a). The finest h = 1.00 mixed- 
element Inflate grid signature is shown with the finest purely-tetrahedral Inflate grid signature in Fig. 22(b). 
The purely-tetrahedral and mixed-element Inflate grid signatures are identical because the grids are identical 
within the core region. The signature are effectively identical for all methods at this close distance, centerline 
location, except for the extrema. 
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The signature at H = 31.8 inches and 0 = 30 is shown in Fig. 23. All four mixed-element grids are shown 
in Fig. 23(a). The signature forward of x = 1.25 remains similar during grid refinement. The signature 
influenced by the model base region x = [1.28, 1.35] changes the most with grid refinement and faceting of 
the sting surface grid aft of x = 1.35 is still evident at this propagation distance. 



(a) Uniformly refined mixed-element grids. 



(b) Various grid methods. 


Figure 23. Delta Wing Body signatures at H = 31.8 and (p — 30. 


The finest h = 1.00 mixed-element Inflate grid signature is shown with the finest purely-tetrahedral 
Inflate grid signature in Fig. 23(b). The BG grid has clustering in the <f> = 30 direction. An example of this 
off-track clustering is shown in Fig. 7(b). The purely-tetrahedral and mixed-element signatures are similar 
with a slight smoothing at x = 1.32 for the Inf-Tet-FUN3D-1.00 and MCAP-Tet-USM3D signatures. 

The signature at H = 31.8 inches and 0 = 60 is shown in Fig. 24. All four mixed-element grids are shown 
in Fig. 24(a). The signature forward of x = 1.25 remains similar during grid refinement. The signature 
influenced by the model base region x = [1-27, 1.35] changes the most with grid refinement and faceting of 
the sting surface grid aft of x = 1.35 is still evident at this propagation distance. 



(a) Uniformly refined mixed-element grids. 



(b) Various grid methods. 


Figure 24. Delta Wing Body signatures at H = 31.8 and cp — 60. 


The finest h = 1.00 mixed-element Inflate grid signature is shown with the finest purely-tetrahedral Inflate 
grid signature in Fig. 24(b). The purely-tetrahedral Inflate, mixed-element Inflate, and MCAP signatures 
are similar with a slight smoothing at x = 1.32 for the purely-tetrahedral Inflate and MCAP signatures. The 
BG grid has clustering in the 4> = 60 direction, which is enabled by the COB improvement to the core grid 
outer surface. A typical vehicle boom carpets is less then 60 degrees, 4 which makes the (j> = 60 extraction 
location uncommon in practice and a robustness test of the grid generation methods. 
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VII. Conclusions 


A review of published methods for specialized near-field sonic boom prediction grid generation methods 
is provided. The influence of uniform grid refinement is studied for a slight (4 degrees) misalignment of the 
grid, where refinement in the propagation direction is more effective than streamwise refinement for reducing 
rounding of the signature. The freestream Mach angle aligned grids exhibited the expected sawtooth wave. 
This study emphasizes the importance of aligning the grid stretching with shocks and expansions for sonic 
boom propagation. Splitting the hexahedral cells of the 3D example into prismatic cells increased the 
rounding of misaligned grids. 

Issues related to the robust extrusion of triangulated cylindrical surfaces is detailed and a description of 
three collar grid extrusion techniques based on tetrahedral core grids is provided. Each method addresses 
the robustness of collar grid extrusion differently. The three methods are applied to models from the First 
AIAA Sonic Boom Workshop to illustrate the utility of these methods. Uniform grid refinement showed grid 
convergence properties of the Inflate method for the SEEB-ALR and Delta Wing Body. These uniformly 
refined grids are required unstructured grids for the workshop. The differences between the signatures 
extracted from the various grid generation methods was equivalent to the magnitude of the experimental 
uncertainty. 
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